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We study multi-photon effects in quantum electrodynamics in lepton photo-production 
on heavy nuclei and nuclear medium at high energies. We are interested in energy, charge 
density and nuclear geometry dependence of the cross sections. We use the impact parame- 
ter representation that allows us to reduce the problem of photo-production to the problem 
of propagation of electric dipoles in the nuclear Coulomb field. In the framework of the 
Glauber model we resum an infinite series of multiphoton amplitudes corresponding to mul- 
tiple rescattering of the electric dipole on the nucleus. We find that unitarity effects arising 
due to multi-photon interactions are small and energy-independent for scattering on a single 
nucleus, whereas in the case of macroscopic nuclear medium they saturate the geometric 
limit of the total cross section. We discuss an analogy between nuclear medium and intense 
laser beams. 

I. INTRODUCTION 

In this article we offer a new perspective on the multi-photon interactions in lepton photo- 
production on a heavy nucleus at high energies. Multi-photon effects arise due to multiple rescat- 
tering of the projectile photon on the protons of the nucleus. Each scattering contributes a factor 
of aZ to the cross section and must be taken into account for heavy nuclei. Traditionally, this 
is accomplished by solving the Dirac equation in the external Coulomb field of the nucleus. An 
explicit assumption of this approach is that the nucleus can be treated as a point-like particle. The 
goal of this article is to solve the photo-production problem in a general case, by explicitly taking 
the nuclear size into account. This allows us to evaluate the finite nuclear size effects in lepton 
photo-production. However, our main interest is to use the developed formalism to investigate 
the photo-production in another extreme case - inside the nuclear medium (e.g. a thin film). We 
will argue that the nonlinear effects are very different in these two opposite cases. The formalism 
that we develop and use in this work can be also applied to the very interesting problem of photo- 
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production on intense laser beams [T]. There the beam size can be tuned so that both limits can 
perhaps be realized experimentally. 

Our approach to the multi- photon effects is based on the Glauber model [5], which allows 
a proper treatment of the nuclear geometry. In Sec. [V] and the Appendix we show that the 
Glauber model calculation of the pair production cross section agrees with the traditional approach 
pioneered by Bethe and Maximon jS], which is based on a solution of the Dirac equation in an 
external Coulomb field (see [31 El El H El O ID] for recent developments). This is an important 
result of this paper. It opens the possibility of a unified approach to the geometric effects in QED. 
In particular, in this paper we study the dependence of the lepton pair photo-production on the 
density of charge sources and nuclear geometry. 

The paper is structured as follows. In Sec. [IT] we write the cross section in a mixed coordinate- 
momentum representation by transforming the scattering amplitude to the transverse coordinate 
space [T2] . The cross section then becomes a convolution of the photon wave function, describing 
splitting of the photon into a lepton-antilepton pair characterized by the transverse vector r, and 
the elastic amplitude iT uz for the scattering of the lepton electric dipole 11 off the nucleus Z at 
some impact parameter B, see Q. The advantage of this representation becomes clear if we recall 
that the trajectory of an ultra-relativistic particle in an external field is a straight line. This implies 
that the dipole representation diagonalizes the scattering matrix. It was shown in [10 that the 
lepton photo-production amplitudes take a rather simple form in the dipole representation. This 
representation was used in [13] for numerical study of the impact parameter dependence of the pair 
production. 



In Sec. Ill we use the Glauber theory for multiple scattering [2] to express the dipole-nucleus 
scattering amplitude iT llz through the dipole-proton one iT llp . This can be done if the nuclear 
protons can be considered as independent scattering centers. For large Z their distribution thus 



follows the Poisson law that leads to a particularly simple expression for the amplitude (10). To 
avoid possible confusion we emphasize, that we consider ultra-relativistic scattering on a "bare" 
nucleus, i.e. a nucleus stripped of all electrons. Such nuclei are provided by the high energy heavy 
ion colliders such as RHIC and LHC. Therefore, the problem of screening of the nuclear Coulomb 
field by electrons - which is an important issue at high energy photon-atom interactions - is not 
relevant here. 



We calculate the scattering amplitude for dipole-proton scattering in Sec. IV. The momentum 
exchanged in the high energy collision is very small. The corresponding impact parameter is 
large. The maximal impact parameter 6 max increases with energy. For a single nucleus, it always 
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holds that 6' max 3> R. The region of large impact parameters is where the dipole cross section 
picks up its logarithmically enhanced energy dependence. Therefore, b' max is the effective radius of 
electromagnetic (EM) interaction of the photon with the nucleus. This is very different from the 
strong interactions where the Yukawa potential drops exponentially at distances ~ l/(2m n ) much 
shorter than the nucleus radius. 

To make connection with the world of high photon densities, we consider photo-production off 
a nuclear medium with transverse extent d S> b' max and the longitudinal extent L <C d. This can be 
thought of as a prototype of an intense laser pulse. Any photon can be converted into an electric 
dipole at the expense of an additional a which, however, impacts only the total rate, but not the 
dynamics of the multiple scatterings. This motivated us to consider the photo-production off the 
nuclear medium in more detail. Also, as has been mentioned, this case is similar in many aspects 
to the well- studied dynamics of the strong interactions. 

At first, in Sec. [V] we study the Born approximation with the leading logarithmic accuracy. 
Born approximation corresponds to scattering on only one proton in a nucleus. For realistic 



nuclei, the photo- production cross section reduces to the Bethe-Heitler formula (31)[l4] with a 
characteristic Z 2 dependence on the proton number. In this case the elastic scattering amplitude 
iT llp is approximately real, see (25), corresponding to the dominance of elastic processes over the 
inelastic ones. This is not difficult to understand: a significant part of the cross section originates 
at distances much larger than the nuclear radius; the nuclear density there is small. On the 
contrary, in the case of the nuclear medium the amplitude iT Up is approximately imaginary, see 



(32), (33c). This implies the dominance of inelastic processes. Formally, this happens because of an 
approximate cylindrical symmetry of the nuclear distribution around a small dipole on one hand, 
and anti-symmetry of the real part of the scattering amplitude Re [ir Hp ] with respect to inversion 



b — > — b, see (12), on the other hand. For the nuclear medium, the photo-production cross section 
is proportional to Z in the Born approximation. 

Multiple rescatterings induce a correction to the Born approximation that generally depend on 



energy ^fs and nuclear charge Z. These effects are studied in Sec. VI Expectedly, the multi-photon 



effects lead to rather different expressions for the resumed cross sections ( 40 ) and ( 43 ) for a single 
nucleus and for the nuclear medium. For a single nucleus the multi-photon processes resumed 
within the Glauber model do not depend on energy. In the Appendix, we demonstrate, using 
the method of [TU] that the Glauber approach reproduces the result of Bethe and Maximon, i.e. 
unitarity corrections beyond the leading logarithmic approximation [3] . In the opposite case of the 
nuclear medium, the cross section grows as a function of energy and Z until it saturates the black 
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disk geometric limit. A problem of propagation of positronium through matter was previously 
discussed by Lyuboshits and Podgoretsky in Ref. [15] who reached similar conclusions. 

A related problem has been recently addressed by C. Muller jTB] who considered scattering of 
an intense laser pulse on a heavy nucleus. He argues - basing on work done with his collaborators 
|17j - that the resumation in the density of the laser beam can be phenomenologically significant 
if one studies the laser-nucleus interaction at the LHC. 

II. PHOTO-PRODUCTION IN THE DIPOLE MODEL 



z 





FIG. 1: Photo-production of a lepton pair. Dashed lines denote photons, solid lines denote fermions. The 
ellipses indicates summation over all proton numbers. Note that the ^-channel photons can hook up to either 
fermion or anti-fermion lines. 

One of the Feynman diagrams for the photo-production process is shown in Fig. [T] In the high 
energy limit s 3> —t, where s and t are the usual Mandelstam variables, the cross section for the 
diagram in Fig. [T] can be written as [18] 

= 2 / Cf £ * = J jf ^kC) J cPBlm(iT lIZ (v,B,s)) . (1) 

This formula is the Fourier transform of the momentum space expression for the scattering am- 
plitude. Dipole size r is a separation of the / and / in the transverse configuration space (i.e. in 
the plane perpendicular to the collision axes); this is the variable that is Fourier conjugated to the 
relative transverse momentum k of the pair. We use the bold face to distinguish the transverse 
two- vectors. Impact parameter B is a transverse vector from the nuclear center to the center-of- 
mass of the 11 pair, see Fig. [2] In ([!]) <£(r, is the square of the photon "wave function" averaged 
over the photon polarizations and summed over lepton helicities. It is given by |12[ ITS] 

c) = 2am2 {^^2 + (1 _ c)2]+ K^rm)} , (2) 
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where m is a lepton mass. £ is the fraction of the photon light-cone momentum taken away by 
the lepton. (iY llz (r, B, s)) is dipole-nucleus elastic scattering amplitude averaged over the proton 
positions (indicated by the (...) symbol) and the dipole-nucleus cross section is, by virtue of the 
optical theorem, 

4f (r, s) = 2 j d 2 Blm {iY llz {i, B, s)) . (3) 

Note, that in the high energy limit, the dipole cross section does not depend on £. Integrating over 
£ we get 

og(s) = {^i 2 H + ^(™)} * lIZ (r,s) . (4) 

Eqs. (|T| and ^ admit the following interpretation. At high energies, the processes of the 7Z 
scattering proceeds in two separated in time stages. The first stage is fluctuation of the photon 
into the lepton pair 7 — ► I + 1. Let the photon, lepton and anti-lepton 4-momenta be q = (<?+, 0, 0), 
k = (£?+, k) and k — q = ((1 — ()q+, (k — q)~, — k), with £ = k + /q + . Here we use the light cone 
momentum notation of the 4- vector = (p + ,p_,p), where p± = p° ±p 3 . The four-scalar product 
is p ■ k = ^{p+k- + p-k+) — p • k. The light-cone time span of this fluctuation is 

A;_ +(q- k)- -q- m 2 +k 2 m 2 +k 2 m 2 + k 2 W 

The largest contribution to the cross section arises, as we will see later, from the longest possible 
Ax+, referred to as the coherence length l c [19]. This corresponds to £ = 1/2 and k = 0, so that 
l c = jptj. Obviously, 1/Z C = Am? /q + is the minimum light-cone energy /_ transfer. 

The second stage of the 7Z process is interaction of the lepton pair with the nucleus. In the 
nucleus rest frame its time span is of the order of the target size 2R, where R is the nucleus radius. 
We observe that l c 3> 2R for sufficiently high g+ for a lepton of any mass m supporting the physical 
interpretation of the two stage process. It now becomes clear why it is convenient to write the 
cross section for the process Fig. [T] in a way in which the photon "wave function" describing the 
splitting process 7 — > I + I is separated from the amplitude Y describing interaction of the 11 
electric dipole with the nucleus. Indeed such an approach proved very useful for the study of the 
strong interactions. 

III. GLAUBER MODEL FOR MULTIPLE SCATTERING 

Let the nucleus quantum state be described by the wave function ipA that depends on positions 
{b a , z a }a=i °f all Z protons, where b a and z a are the transverse and the longitudinal positions of 
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FIG. 2: Geometry of the dipole-nucleus collision in the transverse plane, p, d, and O denote proton, center 
of mass of the dipolc and the nucleus center positions correspondingly, b' = B — b, b < R. 

proton a correspondingly. If the dipole-nucleus scattering amplitude iT uz is known, then averaging 
over the proton positions is performed as 

r z 

(F uz (i,B,s)) = / ]Jd\ a dz a \^i,^2,^,...)\ 2 T llz (v,B-b 1 ,z 1 ,B-b 2 ,z 2 ,...,s). (6) 

J a=l 

The scattering amplitude is simply related to the scattering matrix element S as T(s, B) = 1 — 
S(s,B)- The later can in turn be represented in terms of the phase shift x so that in our case 

T lIz (i, B - bi, zi, B - b 2 , z 2 , . . . , s) = 1 - exp{-i/ z (r, B - bi, z Xl B - b 2 , z 2 , . . . , s)} . (7) 

At high energies, interaction of the dipole with different protons in the nucleus is independent 
inasmuch as the protons do not overlap in the longitudinal direction. This assumption is tanta- 
mount to taking into account only two-body interactions, while neglecting the many-body ones [2]. 
In this approximation the phase shift x i n the dipole-nucleus interaction is just a sum of the 
phase shifts x llp m the dipole-proton interactions i.e. 

T iIz (t, B - bi, -si, B - b 2 , z 2 , . . . , s) = 1 - exp{-i ^ X lIp (r, s,B- bj} . (8) 

a 

The problem of calculating the scattering amplitude of the dipole on a system of Z protons thus 
reduces to the problem of calculating of the scattering amplitude of the dipole on a single proton 
T llp . In this approximation, (e~ tx ) = e~ llyX \ i.e. correlations between nucleons in the impact 
parameter space are neglected. 
We can re- write © and (|7j) as 

<r^(r, b, *)) = l - (lit 1 - r " P fe s ' B - b a )]> • (9) 

a 

For a heavy nucleus with large Z we obtain 

(T llz (r, B, s)) « 1 - ex P {-Z(r^(r, s,B — bj)} . (10) 
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This approximation is justified as long as a -C 1 and aZ ~ 1. Indeed, expression in the exponent 



of (10) can also be obtained by expanding 



flz - Zln(l - r" p ) « -Z[T lTp - -(T lTp ) 2 + ..]~aZ + 0{a 2 Z) 
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IV. DIPOLE NUCLEUS SCATTERING AMPLITUDE 

Now we turn to the calculation of the dipole-proton amplitude, see Fig. [3] At high energies 
the t-channel Coulomb photons are almost real: I 2 ~ —I 2 . This is the well-known Weiszacker- 
Williams approximation that we use throughout the paper. It allows us to determine the leading 
logarithmic contribution to the scattering amplitude. In this approximation, the real part of the 
elastic dipole-proton scattering amplitude reads: 



Re[ir" p (r,b',s) 



a 


fcPl r 


«J 


' P . 


2a 


r°° di 






Jo T 



e <(b'+|r).l _ e i(b'-ir).r 
[Jodb' + irlO-Jodb'-irlO]. (11) 



Taking the integral we arrive at 



lb' - irl 

Re [ir" p (r, b', s)] = 2a In ^ f=[ . (12) 



lb' + |r| 



The nuclear density p(b, z) is normalized such that 

J d 2 bdzp(b,z) = Z. (13) 

Therefore, the average of the real part of the amplitude over the proton position reads 

1 f f z ° lb' — -rl 

<Re[ir"P(r,B,s)]> = - d 2 b d zp2aln l —^\, (14) 
Z J J-zq |b + |l| 

where zq = 2\/ R 2 — b 2 . Since p ~ const, we get 

(Re[ir»f(r,B, S )]) = ^ / d 2 b P 2^W^b 2 In , (15) 

z J |b + p| 

with b' = B - b, see Fig. [§ 

The imaginary part of the dipole-proton elastic scattering amplitude is depicted in Fig. |3 
The corresponding analytical expression reads: 

y 2 rd 2 i rd 2 r 



Im^(r,b', S )] = ^/^/^ (e'(b'+ir)-l_ e i(b'-§r)-l) (^(b'+Ir).!' _ e -i(b'-|r).l' 



(16) 
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Averaging over the proton positions we derive 
(Im[ir»f(r,B, S )]) = iy d 2 6 P 2Vi? 2 -6 2 ^ 

21/ 



The Weiszacker- Williams approximation breaks down for very small momenta of the t-channel 
photons. To determine the cutoff transverse momentum Z m j n at which this approximation still holds, 
we use the Ward identity applied to the lower vertex of Fig. |3ja). It gives p- 1 = p+l- +p-l+ = 0, 
where p = (M 2 /p_,p_,0) with M being the proton mass. It follows that l + = —p + l^/p-. Thus, 
—I 2 = + 1 2 w l 2 if l 2 3> (p+/p-)Z 2 . Since the minimum value of /_ is 4m 2 /q + (see ^ and 

the following text) we find 



p+ (2m) 4 _ M 2 (2m) 4 _ M(2m) 2 

\J p- q+ V piq+ s 

where s = (p + q) 2 as p~q+. ludn corresponds to the maximum possible impact parameter 6' max ~ 
1/^min- As we will see in the next section, the energy dependence of the dipole-nucleus amplitude 
arises from the logarithmic dependence on these cutoff scales (this is why a precise value of the 
proportionality coefficient between b' max and l/7 m in is not very important). 

V. BETHE-HEITLER LIMIT 

Before we consider a general case it is instructive to discuss the scattering process in the Born 
approximation. We would like to (i) confirm that we reproduce the result of Bethe and Heitler 
|14j and (ii) determine the range of impact parameters B that gives the largest contribution to the 
cross section. 

The Born approximation corresponds to taking into account only the leading terms in the 



coupling a. Expanding the exponent in (10) we have 



(T iIz (t_,B,s)) = Z(T iI p(t,B,s)} -^Z 2 (T iIp (t,B,s)) 2 + ... . (19) 

The real part of iT llp stems from a diagram with a single photon exchange, while its imaginary 
part arises from the double photon exchange (in the high energy limit these diagrams are purely 
real and imaginary correspondingly), see Fig. Therefore, at the lowest order in a, the imaginary 



part of iT llp corresponds to the two photon exchange Fig. pfb) that gives - according to (19) 



contribution of order a 2 Z ~ a to the cross section rt3J). On the other hand, the real part of (iT llp ) 



\2 



9 



iT il P = | I 

v-L- 



Re\iT lf P} ^Im[^^ ^^ P] 
(a) (b) 

FIG. 3: The first two terms in the expansion of T Up in a. 

corresponds (at the leading order in a) to the product of the two diagrams Fig. [3^a) with single 
photon exchanges; it is of the order a 2 Z 2 ~ 1 *. Thus, in the Born approximation 



ag B (r, S ) = 2 j d 2 B j^Z 2 (Re [iT" p (r,B, s)}) 2 + Z (Im [iT llp {i, B, s)])| . 



(20) 



The reason for keeping a sub-leading term in (20) (the term in the integrand proportional to Z) 



is that in the nuclear medium the real part of the dipole-proton amplitude vanishes, see Eq. (|32j) 
and then the main contribution stems from the imaginary part of that amplitude. 

A. Single nucleus 

Values of the dipole-proton impact parameter b' can range up to its maximal value 6' max , which 
is determined as the inverse minimal momentum transfer, i.e. b f max ~ 

lAmin, where Z min is given by 



(18). For electron production at the collision energy of ^fs = 10 GeV we estimate b' max ~ 20 nm 
whereas the radius of the Gold nucleus is 7 fm. Thus for realistic nuclei & max 3> R. (In the case of 
r photo-production this approximation holds for energies larger than about yfs = 30 GeV). 

We will demonstrate shortly, that the large logarithmic contribution to the dipole cross section 



(20) comes from the region b' 3> b and b' S> r/2. Because by definition b < R, this is equivalent 
to B ~ b' S> max{i?, A/2}, where A = 1/m is the Compton wavelength of a lepton. Here we used 
the fact that the characteristic value of the dipole size is r ~ 1 jm since it corresponds to a sharp 
maximum of the integrand of Q. Note, that for e~e + dipole R <C A <C 6' max , while for /i and r 
dipoles A < R < b' max . 



* This is in sharp contrast with QCD where the single gluon exchange amplitude is identically zero because the non- 
Abelian symmetry generators are traceless. Thus, there is no elastic scattering at this order. As a consequence 
the leading contribution in QCD is of the order of atomic weight A. 
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Using b' B 3> b, see Fig. [2] we can write ( 15 ) as 

IB - k| 1 



(Re[ir e e+ P(r,B,g)]) « 2 a In ^ - f ~ ^ d 2 b p2^ R 2 - b 2 = 2 a In 



|B - ±r| 
IB + |l| 



(21) 



To obtain the cross section we need to integrate over the impact parameter B. Expand the logarithm 

IB ~h\ B • r 



in at £ > r/2 as 



In 



+ 0(r 3 /5 3 ) 



B + ir| B 2 



(22) 



we derive 



(Re[ir e e + P(r,B,s)]) «2a 



Br 



B 



2 • 



(23) 



Turning to the imaginary part we can write Eq. (17) as follows 

„2 



(Im[ir e+e ^(r,B, S )]> = ^ 



x e 



2;/ 



(B+|r).l _ J(B-|r).n / -i(B+lr).l' _ p -i(B-ir).l' 



4 a 2 In 2 



IB 



|B + |r| 



2tf 



(B-r) s 



(24a) 
(24b) 



To the leading order in a we can neglect the small absorption effect given by (24b) and write for 
the scattering amplitude 

Br 



(r"P(r,B,s)> ^i2a 



B 



(25) 



Consequently, in the high energy and Born approximations the dipole-nucleus cross section becomes 



a l t l l B (v,s) = Z' I d A B4a 



(B • r) 2 



(26) 



Denote 1Z = max{i?, A/2}. Integration over B in the interval max{i?, r/2} < B < 6 max gives 



cT l t l l B (i,s)=4na 2 r 2 Z 2 ln 



2 2y2i "max 



n 



(27) 



where the logarithmic accuracy allowed us replacement r = 1/m in the argument of the logarithm. 
The total photo-production cross section off a nucleus at the order a 3 in the perturbation theory 
reads: 

d 2 r <- 1 



7 z 

tot,£ 



T 

2^ 



1 



/>oo 

V4a 3 Z 2 m 2 / drr 
Jo 



'-K\{rmi) + Kl(rmi) 



2 1 "max 

r In 



n ' 



(28) 
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where the sum runs over all lepton species. Taking integrals over the dipole sizes finally produces 

a ]Z ( 8 ) = Y4a 3 Z 2 m* (-— A + In ^ . (29) 

tot.BW l_j 1 \3 3mf 3mfJ K v ; 



In particular, for electrons 



yZ^e+e-X, \ 28 a 3 Z 2 s . . 

cr' R (s) = ?^ln . (30) 

tot,B V ) 9 o 2Mm e V 7 



In a frame where nuclear proton moves with velocity v = (p+—p-)/(p++p-) the minus component 
of its four-momentum is p_ = My/(1 — v)/(l + v). Therefore, s = 2Mwy / (1 — v)/(l + v), where 



the incoming photon frequency is oj = q+/2. Eq. (30) now reads 

- t 7 otB £+e ^(*) = ^ "Tin f ^.Z 1 ^! • (31) 
t *- 8 w 9 m? e \ m e V 1 + 

This is the high energy limit of the formula derived by Bethe and Heitler |14J . In the Appendix 
we discuss the energy-independent correction to this equation. 



B. Nuclear medium 

Consider a nuclear matter having an arbitrary large number of protons Z. We assume that the 
inter-proton distance is much smaller than the Compton wavelength of electron, i.e. p -1 / 3 <C A. 
This allows treating the nuclear matter as the continuous medium of density p. Let the medium 
transverse size be d > b' mSLX and longitudinal size L <C l c - In the medium rest frame l c ~ b' maJC 
Therefore, the nuclear medium has a form of a film with d 3> L. What does the lepton photo- 
production look like in this case? The key observation in this case is that the dipole-proton 
interaction range is much smaller compared to the transverse size of the medium. Therefore, 



b' <C B and then it follows from ( 15 ) that the real part of the elastic dipole-proton scattering 



amplitude vanishes. Indeed, changing the integration variable b — > b' in (15) and taking the limit 
of small b' we obtain 

2a f ,,,, ,„ , , s _ , |b' - |r| 



(Re[ir^(r,B, S )]) = ^jr J d 2 b'p(B - W) L In ^ 



b' + ^r 



2a p{B)L /d 2 6'ln^-M = 0. (32) 
J |b + p| 



z 

The last equation follows because the integrand is an odd function under the reflection b' — ► — b . 
The first non- vanishing contribution to the real part of the amplitude is proportional to the nuclear 



' The values of b' m!ix at yfs = 10 GeV for e, /i and r are 20 nm, 0.4 pm, and 2 fm correspondingly. 
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density gradient, which has the largest value at the boundary. These small effects can be taken 
into account if the nuclear density distribution is known. 

The leading contribution arises from the imaginary part of the amplitude that reads 



<Im[^(r,B, S )]) « ipL I d 2 b' ^ 



i(b 



: /f(-1 

2 



z pL4a' J ^ 

— pL4a 2 (2vr) -r 2 ln — 
Z 4 rZ 



(33b) 
(33c) 



where in the last line we kept only the logarithmically enhanced term. We see that in the case of 
the nuclear medium, the entire elastic dipole-proton amplitude is approximately imaginary and is 



given by Eq. (33c). Moreover, unlike (24b) it increases logarithmically with energy. Substituting 



(32) and (33c) into (20) we obtain in the Born approximation 

2 



a l t ^ B (r_,s)=4nZa 2 r 2 ln 



(34) 



Note, that cross section (34) is Z times smaller then the one given by (27). 



VI. UNITARITY EFFECTS 



The unitarity relation applied to the elastic scattering amplitude at a given impact parameter 
reads (we suppress the argument (r, b, s) of all functions) 



2Im(ir" z ) 



+ G 



uz 



where G llz is the inelastic scattering amplitude. Using (1 1 Oh we can solve (|35|) as 



G 



UZ e -z<im(ir"p)> 



'tot 

jTz 



uz 

°el 



(35) 



(36) 



It follows that the total, inelastic and elastic cross sections are given by 

oE = 2 J d 2 B {l - cos[Z(Re (i^))] e ~^^ lTp ))^ (37a) 

J d 2 B[l-e- z ^ lXp ))) (37b) 

J d 2 B{l-2cos[Z(Re(ir"P))]e~ z < Im ( jr " P »+ e ~ z < Im ( ir " P »} . (37c) 

When Z(T ll P) ~ 1 deviation fr om the Born approximation becomes large and the unitarity 
effects set in. We proceed to analyze the unitarity corrections in two extreme cases. 
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A. Single nucleus 



In the case of photon-nucleus scattering, we derive using (25) and integrating over the directions 
of B: 



B • r \ [Bmsx 

% rr.-.j = 2 / d 2 B ( 1 - cos{2Z«^}j = 4vr J dB B [1 - J {2Zar/B)} . (38) 



Introducing the dimensionless variable x = 2Zar / (2B max ) we can express this integral in terms of 
the generalized hypergeometric function: 



4! (E, «) = (2Z " r)V {x 2 2 F 3 [(1, 1), (2, 3, 3), -x 2 ] + 4(2 - 2 7 - In* 2 )} . 



(39) 



where 7 is Euler's constant. The total photo-production cross section of light leptons is obtained 
using Q . Actually, it is convenient to directly plug ( 38 ) into Q and first integrate over r and then 
over B. The result can be expressed in terms of another dimensionless variable y = aZ/(B max m) 
as: 



g 3 Z 2 
9m 2 y 3 



2y-28y 3 ln(2y) 



+{-3^1 + y 2 + 24y 2 v / l + 2/ 2 + cosh[3 ln(y + \A + y 2 )]} ln(y + ^1 + y 2 



(40) 



To find the high energy asymptotic we expand (40) at small y and find 

28a 3 Z 2 



1 



41 12 



9m 2 



■[In 1 1 w z + 

L 2y 42 35 y 



(41) 



In the leading logarithmic approximation that we use in this paper, only the leading logarithm in 



the square brackets of (41 ) can be guaranteed. Indeed, it reproduces the Bethe-Heitler formula (30 ). 



However, Eq. (41) does allow us to conclude that corrections to the leading logarithm must not 



increase with energy. Therefore, in the case of lepton pair photo-production on a single nucleus 
the unitarity corrections do not grow with energy at the leading order in a <C 1 and aZ ~ 1. 
Actually, the Glauber model that we use in this paper also allows us to reproduce the subleading 



(i.e. energy-independent) corrections to Eq. (30) that were first obtained by Bethe and Maximon 



[3]. This result can be easily derived using the formalism developed by Ivanov and Melnikov in [10 . 
The derivation is outlined in the Appendix. In the framework of our approach it is straightforward 
to calculate corrections of order a 2 Z to the Bethe-Maximon formula by taking into account the 



imaginary part of the amplitude (24b). This problem will be addressed elsewhere. 
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B. Nuclear medium 



Conclusions of the previous subsection are dramatically reversed in the case of scattering off a 

into 



nuclear medium. For the total dipole cross section we obtain after substituting (32) and (33c 



(37) 



o-f f(r,s) = 2ttgP(1 -exp{-pL27raVln 



2 2 ■ 



rl. 



-}) 



Similar results were obtained for an infinite medium in [15]. Substituting into Q we derive 

d 2 r (2 



°tot 0) = 2a m d 



-K 1 (rm) + Kq (rm) 



1 



-pL2ira 2 r 2 ln 



(42) 



(43) 



2vr |3 

At large Z the multi-photon effects start to play an important role. 

It is instructive to study the unitarity corrections by increasing the proton number Z while 
keeping the medium size given by L and d fixed. This corresponds to an increase in the proton 
number density. We are interested to know at what Z's the multi-photon effects become observable. 
In Fig.|4] we show the ratio of the total photo-production cross calculation to its Born approximation 
for = 10 GeV and L = 1 nm, d = 100 nm. We observe that for electrons the deviation from the 
linear Born regime starts at about Z ~ 10 12 charges. At Z ~ 10 13 the cross section is essentially 
black (i.e. Z- independent). 




FIG. 4: Ratio of the lepton photo-production cross section and its Born approximation as a function of the 
"nucleus" charge Z. y/s — 10 GeV, L=lnm, and d=100 nm. 



VII. DISCUSSION 

In this article we discussed an effect of multiple scattering on the high energy lepton photo- 
production off a heavy nucleus and nuclear medium. The heavy nucleus serves as a source of strong 
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Coulomb field. Resummation of multiple scattering of a lepton dipole in the nucleus amounts to 
resummation of diagrams with an arbitrary number of photons attached to the lepton loop. In the 
heavy nucleus case we selected only diagrams that are enhanced by the power of Z and resummed 
them using the Glauber model. We found that dominance of elastic scattering leads to energy- 
independent unitarity corrections. These corrections do not saturate the geometric limit of the 
cross sections. The unitarity corrections that grow with energy can appear as a result of multiple 
lepton pair production at very high energies [20 . The dipole model that we used in this paper is 
a convenient tool for taking into account such high energy evolution effects. We plan to address 
this problem elsewhere. 

Lepton photo-production off a nuclear medium is different in many aspects. Here the processes 



are predominantly inelastic. At very large Z the cross sections (37) saturate at the geometric limit 
due to multi-photon interactions with the medium. An additional contribution to the saturation 
stems from the multiple electron-positron pair production |20|, [21] . This leads to growth of the 
imaginary part of iT llp as a power of energy s A with A ~ (ll/32)7ra 2 [2T]. However, since A is 
very small this effect is mostly of academic interest. Parametrically, lepton photo-production off 
the nuclear medium is very similar to the hadron photo-production in strong interactions. The real 
part of iT llp vanishes identically in that case because it is proportional to the vanishing trace of the 
Gell-Mann matrices. In QCD, the quark-anti-quark color dipole interacts with the heavy nucleus 
by means of multiple exchange of virtual gluons. Additionally, multiple gluon emission contributes 
to the power-law dependence of the cross section. The corresponding power is not small leading to 
important observable effects. Unitarity correction in QCD grow very large at high energy leading 
to saturation of the scattering amplitudes |22[ [23], [23] . In lepton photo- production QCD evolution 
will arise at high orders of perturbation theory. 

In the case of the nuclear medium, weakness of the EM coupling can be compensated by large 
number density of charges. In a more realistic model one should also take into account screening 
of the Coulomb fields by electrons leading to the emergence of the nuclear form-factor. However, 
our emphasis in this article is the study of electro-magnetic interactions in a high density system. 
A high density of charges can be realized also in ultra-fast pulse lasers that emit pulses containing 
as many as N ~ 10 20 photons and having dimensions L < d. In this case the inelastic in-medium 
effects start to play a significant role in large range of collision energies. Although fluctuation 
of laser photons into electric dipoles is suppressed by the power of a, the created Coulomb field 
is enhanced by large effective charge Z. We are going to address this intriguing problem in a 
forthcoming publication. 
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APPENDIX A: DERIVATION OF THE BETHE-MAXIMON FORMULA 



Eq. (38 ) is written in the leading logarithmic approximation that does not allow one to reproduce 



the subleading energy-independent terms in the total cross section. However, the Glauber model 
allows one to calculate the total cross section with better accuracy. To this end, one has to abandon 



the approximation of (22) and use formula (21) for the real part of the amplitude in (37). We have 

,2 



°tot( s ) 



am 



TT 



— \-K 2 (mr)+K 2 (mr)}2 [ d 2 B{l - cos \2aZ\n | B l/2 \ 
2ir L 3 1V ' 0V n J \ L |B +r/2| 



(Al) 



Introducing a new integration variable R = B + r/2 we re- write ( Al ) as 



-fZi (Mir I d 2 r T 2 



Consider 



TT 



\-K?(mr) + K. 
2tt l 3 11 ' ' 



2 (mr)] 2 J 



d 2 R\l 



1 /IR 



R 



2iaZ 



1 /|R-r| 
2V R 



-2iaZ 



Aa 



7 z 

tot 



'tot 



; tot,B 



(A2) 



(A3) 



i.e. the contribution due to multiple rescatterings. The term corresponding to the Born approxi- 



mation is obtained by expanding the expression in the curly brackets in (A2) at small aZ. One 
then arrives at the formula derived by Ivanov and Melnikov |10j . It reads: 
- 2 '-d 2 r r 2 



tot 



am 



IT 



^ [-Kf(mr) + Kl{mr)} 



d 2 R< 1 



l /IR 



2 V R 

Taking integrals over R and r yields: 



2iaZ 



1 / IR - rl 
2 



R 



-2iaZ 



4(aZ) z In 



2,_2 lB-E| 



R 



Aa. 



28 ^v^\ 
-a Z 



[1 - iaZ) + V(l + iaZ) - 2^(1)] , 



(A4) 



(A5) 



9m 2 " 2' 

where ip{z) is the logarithmic derivative of the Gamma function. As expected, the unitarity 



corrections are independent of energy in agreement with the results of Sec. VI A 
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To calculate the energy-independent term of the Born approximation (the last term in the curly 



brackets of (A4)) one needs to take into account the fact that the smallest longitudinal momentum 
transfer depends on £ and k, see ([5]). In the logarithmic approximation, we maximized ([5| with 
respect to both £ and k. Now it must be retained which means that (i) the regularization of the 
logarithmically divergent integrals must be done in momentum space and (ii) the £ integral (see 
Q) must be done after the regularization since the cutoff depends on it. We refer the interested 
reader to the Ref. [10] for details. The result is 

28 o?Z 2 



1 z 



; (s) 



In 



109 \ 
42~ J ' 



(A6) 



tot,sv-/ c, ,„•> \ ' " Mm e 
The second term in the brackets is the subleading correction we were looking for. The final result 



is the sum of ( A5 ) and ( A6 ) in agreement with |3 . 
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